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ABSTRACT 



Context. Several clusters of red supergiants have been discovered in a small region of the Milky Way close to the base of the Scutum- 
Crux Arm and the tip of the Long Bar. Population synthesis models indicate that they must be very massive to harbour so many 
supergiants. Amongst these clusters, Stephenson 2, with a core grouping of 26 red supergiants, is a strong candidate to be the most 
massive cluster in the Galaxy. 

Aims. Stephenson 2 is located close to a region where a strong over-density of red supergiants had been found. We explore the actual 
cluster size and its possible connection to this over-density. 

Methods. Taking advantage of Virtual Observatory tools, we have performed a cross-match between the DENIS, USNO-B1 and 
2MASS catalogues to identify candidate obscured luminous red stars around Stephenson 2, and in a control nearby region. More 
than 600 infrared bright stars fulfill the colour criteria, with the vast majority having a counterpart in the / band and > 400 being 
sufficiently bright in / to allow observation with a 4-m class telescope. We have observed a subsample of ~250 stars, using the multi- 
object, wide-field, fibre spectrograph AF2 on the WHT telescope in La Palma, obtaining intermediate-resolution spectroscopy in the 
7500-9000A range. We derive spectral types and luminosity classes for all these objects and measure their radial velocities. 
Results. Our targets turn out to be G and K supergiants, late (>M4) M giants, and M-type bright giants (luminosity class II) and 
supergiants. We find ~35 red supergiants with radial velocities similar to Stephenson 2 members, spread over the two areas surveyed. 
In addition, we find ~40 red supergiants with radial velocities incompatible in principle with a physical association. 
Conclusions. Our results show that Stephenson 2 is not an isolated cluster, but part of a huge structure likely containing hundreds 
of red supergiants, with radial velocities compatible with the terminal velocity at this Galactic longitude (and a distance ~6 kpc). In 
addition, we find evidence of several populations of massive stars at different distances along this line of sight. 
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1. Introduction 

In the past few years, a flurry of discoveries has revealed several 
clusters of red supergiants (RSGs) located in a smal l region of 
the Galactic plane, between I = 24° and / = 29° (e.g.,|F iger et al 



2006; iDaviesetalJ 12001 IClarket al.l [2009ah iNegueruel a et al 



2011). These clusters are so heavily obscured that, until now, 



the only members observed are the RSGs. Population synthesis 
models suggest that the clusters must conta in very large stella r 
populations to harbour so many RSGs (e.g jDavies et alj|2007l) . 
Recent simulations indica te that, at typical ages r = 10-20 Myr 
dDavies et aDl2007L [20081). one should expect ~10 4 M Q in stars 
for each 3-5 RSGs Clark et alj|2009bh . 

RSGC2 = Stephenson 2 (Ste 2; / = 26°2, b = -0°1) is 
the least obscured and apparentl y most mass i ve of the red su- 
pergiant clusters. Discovered by Ste phenson! dl990[). the clus - 
ter was foun d to ha ve a clump of RSGs ([Ortolani et al. 2002). 
iDavies et al.l (I2007L henceforth D07) obtained A"-band spec- 



troscopy of a large number of bright sources within 7' of 
the nominal cluster centre, finding more than 25 RSGs that 
shared similar radial velocities. After defining an average ra- 
dial velocity for the cluster and eliminating outliers, they found 
vlsr = +109.3 + 0.7 kms -1 , with the uncertainty represent- 
ing Poisson statistics for 26 likely members with measured val- 
ues (|D07). Stars were considered members if their vlsr was 
within +10 kins" 1 of the average value. On the other hand, 
iDeguchi et al.l (2010a) obtain vlsr= +96.5 kms 1 from mea- 
surements of SiO masers associated with four cluster members. 
They attribute the difference to systematic effects. 

The age of Ste 2 has been estimated at t= 17 + 3 Myr 



Ortola niet al.l 



D07 



estimate, 



* Based on observations collected at the William Herschel Telescope 
(La Palma, Spain) 



(D07), consistent with a value ~20 Myr derived by 
(120021) . Assuming that the 26 RSGs are members, 
using population synthesis models, that the und erlying cluster 
popul ation must be M c \ 5 x 1 4 Mq . Converselv. lDeguchi et al.l 
(l2010al) claim that the spatial distribution of RSGs argues for 
the existence of two clusters, which they call Ste 2 and Ste 2 
S W, though their radial velocities do not seem to be significantly 
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different. In summary, the actual spatial extent, and hence total 
mass, of Ste 2 is still poorly determined. 

The line of sight to the RSG clusters (RSGCs) passes through 
several Galactic arms and reaches an area of very high ob- 
scuration at a distance comparable to those of the clusters 
( Negu eruela et al ] l2010h . Stellar densities are extremely high at 
moderate K magnitu des, a nd membership of individual stars 
is difficult to assess. ID07I found a large population of field 
stars with Ks magnitudes comparable to the supergiants. Indeed, 
down to their limiting magnitude Ks = 8, there are many more 
field stars than RSG members in the r = T circle analysed 
(thou gh RSGs are more numerous than field stars for Ks < 6.2). 
|D07| identify these field stars as predominantly foreground gi- 
ants, though several putative RSGs with v ra( j suggestive of non- 
membership were found. 

In a recent paper (Negu eruela et al.l 1201 it from now on, 
Paper I), we studied the spatial extent of another cluster of 
red supergiants, RSGC3, w ith estimated t= 16-20 Myr and an 
inferr ed M c i=2-4xl0 4 M o (IClark et al.ll2009at lAlexander et alJ 
2009), finding that it is part of an association, which in- 
cludes smaller clusters of red supergiants. This association 
contains <: 50 RSGs and hence a total mass > 1O 5 M . 
This size has been confirmed (and likely surpassed) by the 
recent discovery of another cluster of RSGs in its vicinity 
(Gonzalez-Fernandez & Negueruela 2012). 

The connection of the RSGC3 association to Ste 2 and the 
other RSG clusters is unclear. Though all the clusters are lo- 
cated at a similar distance and have similar ages, they span 
~500 pc in projection. There have been suggestions that this 
Scutum Complex represents a giant star formation region trig- 
gered by the dynamical excitation of the Galactic Bar, whose 
tip is believed to inters ect the Scutum-Crux Arm close to this 
region (see discussion in lD07tlGarzon et al ]|1997l) . The possibil- 
ity of a hierarchical star-forming structure extending over such 
a long distance seems unlikely as it would be larger than the 
largest cluster com plexes seen in exte rnal galaxies (for instance, 
~ 240 pc in M5 1 ; Bastia net al.l2005l) . even t hough less cohe rent 
structures (aggregates in the terminology of lEfremovll 1 9941) are 
seen spanning comparable distances. 

All the RSG clusters have similar ages, in the 10-20 Myr 
range, but this do es not guarantee an association. RSGC1, with 
an age -12 Myr dDavies et alJ l2008). could be more closely al- 
lied to its close neighbour, the Quarte t cluster, with an age be- 
tween 3-8 Myr dMessineo et al.ll2009l) . while the RSGC3 asso- 
ciation lies very close in the sky to the W43 star-forming com- 
plex, also at d ~ 6 kpc, whi ch covers the I = 29°.5 - 3 1°5 range 
(iNguven Luong et alJl201ll) . 

The distribution of Galactic Hn regions shows a 
marked maximum towards I ~ 30°-31° (Ba nia et alJ [2010: 
lAnderson et al.ll201 lb . where tens of sources display radial ve- 
locities in the vlsr ~ 90 - 105 kms -1 interval. This is generally 
interpreted as the tangent point of the Scutum-Crux arm. There 
is a secondary maximum towards / ~ 24°, where many sources 
have vlsr ~ 110 km s -1 , similar to the systemic velo city of 
RSGC1 (/= 25°.3; v LSR = 123 km s'^ lDavies et alJl2008h . There 
is a local minimum towards / ~ 26°, th ough the number of H n 
regions in that direction is still high (lAnderson et al.l 1201 ll) . 
Gala ctic molecular clouds s eem to follow a similar distribution 
(e.g jRafhborne et al.ll2009h . 

In this paper, we study the distribution of luminous red stars 
in the vicinity of Ste 2. Though the cluster core is clearly defined 
by a strong concentration of bright infrared sources, its actual 
extent is poorly determined. The spatial distribution of known 
members is very elongated, and the possibility that the cluster 
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Fig. 1. The spectrum of one of our targets (J18414143-0531378) 
observed with AF2 (top), at R = 4 100 and with ISIS at R » 
16 000. At the AF2 resolution, many of the lines seen in the 
ISIS spectrum are blended into broad features. An interesting 
example is the 8514A blend (mainly Fei) that appears as a sin- 
gle feature at the AF2 resolution. This feature is very strong in 
luminous supergiants (see Fig. Et- 



is surrounded by a halo has not been studied. Likewise, the fact 
that RSGC3 is surrounded by an extended association opens the 
possibility that Ste 2 is not isolated eit her. 

Following the strategy outlined in iPaper J. we have used a 
multi-object spectrograph with high multiplexing capability to 
obtain intermediate-resolution spectra of a large sample of bright 
infrared sources in the spectral region around the Ca n triplet. In 
Section [2] we discuss the photometric criteria used for selection 
of candidates and describe the observations of a sample of these 
candidates in the region around Ste 2 and a second sample in 
a control area about 1° away. In Section [3] we summarise the 
analysis methodolog y. We ha ve then used the criteria and cali- 
brations developed in lPaper II to assign spectral types and calcu- 
late colour excesses. We also measure radial velocities for all our 
sources. The results of the analysis are presented in Section [4] 
Finally, we discuss their implications in Section and wrap up 
with the conclusions. 



2. Target selection and observations 

In IPaper B we showed that a careful application of photometric 
criteria over existing catalogues could select obscured red lumi- 
nous stars with high accuracy. Unfortunately, the photometric 
properties of RSGs do not allow a clear separation from M-type 
giants, implying the need for spectroscopic classifications. 

2. 1 . Selection of candidates 

Three catalogues have been used to search for candidates. 
We have started with the 2MASS Point Source Catalogue 
dSkrutskie et all l2Q06h . 2MASS contains photometric data for 
471 million sources in the J (1.25/im), H (1.65/im), and Ks 
(2.16^m) near-infrared bands. Typical photometric and astro- 
metric uncertainties of 2MASS are less than 0.03 mag and 
100 mas, respectively. Three selection condition s were applied 
on these data, based on the criteria discussed in Paperlt 

1. We selected stars bright in the Ks band: Ks < 7.0. 

2. We selected reddened stars by requiring a colour redder than 
that of any star: (/ - K s ) > 1.3. 



2 



I. Negueruela et al.: Red supergiants around the obscured open cluster Stephenson 2 



3. We selected stars with infrared reddening-free parameter, 
Qir = (J - H) - 1.8 X (H - K s ), typical of RSGs: 0.1 < 
Qir < 0.4. Though some RSGs have £>ir outside this range 
(cf. lMessineo et al.ll20T2h . this cut selects most RSGs, while 
leaving out low-luminosity red stars. 

The vast majority of our targets present a good quality flag 
(Qflg = A, B, C, or D) in all three 2MASS bands. A few ob- 
jects with other quality flags were included in the sample and 
assigned low priority, after we checked out that they were not 
badly blended. In the end, only two such objects were observed, 
J 1 8402969-05463590 and J 1 8403226-06 1 1584. 

The 2MASS objects selected were subsequently cross- 
matched with two others catalogues, DENiSQand USNO-B.l 
dMonet et al.ll2003l) . DENIS provides photometry in the Gunn ; 
band (0.8/vm), while USNO-B.l gives magnitudes in a photo- 
graphic near-infrared "7" band. We applied two magnitude cri- 
teria to the objects cross-matched: i) we rejected all stars bright 
in I (I < 10), so as to leave out unreddened stars; and ii) we re- 
quired / < 17, to select only stars observable with current instru- 
mentation. We found more than 600 candidates within 1° of the 
centre of Ste 2 in the 2M ASS/DENIS cross-match, and a slightly 
smaller number in the 2MASS/USNO cross-match. Around 30% 
of the USNO sources do not have counterparts in DENIS, result- 
ing in slightly above 800 unique candidates. These are not all 
the candidate red luminous stars in the area, but only those with 
a counterpart in the far red with 10 < I < 17. 

For the AF2 observations presented in this paper, we had to 
restrict the target list to objects brighter than I = 15.1, to make 
sure that the star would give a usable spectrum with a 4 m class 
telescope. Using this new criterion, the main target area (defined 
below and covering ~ 0.8 deg 2 ) contains 226 2MASS/DENIS 
cross-matches, and 202 2MASS/USNO cross-matches, with 246 
unique objects. 

To apply the colour-magnitude criteria and to cross-match 
the catalogues, we took advantage of the tools offered by the 
Virtual Observatory (VO), an international project designed 
to provide the astronomical community with the data access 
and the research tools necessary to enable the scientific ex- 
ploration of the digital, multi-wavelength universe resident in 
the astronomical da ta archives. In particular, we used Aladin 
dBonnarel et alj|2000T) . a VO-compliant interactive sky atlas that 
allows users to visualise and analyse digitised astronomical im- 
ages, astronomical catalogues, and databases available from the 
VO services. We also used TOPCAT0, another VO-tool, to purge 
duplicated data. TOPCAT is an interactive graphical viewer and 
editor for tabular data that allows the user to examine, analyse, 
combine and edit astronomical tables. 

2.2. Observations 

Observations were carried out with the AutoFib2+WYFFOS 
(AF2) multi-object, wide-field, fibre spectrograph mounted on 
the Prime Focus of the 4.2 m William Herschel Telescope 
(WHT), in La Palma, Spain. The observations were taken on the 
nights of 2009, June 5th (in service mode) and June 6th-7th (in 
visitor mode). AF2 uses 150 science fibres with diameter l'/6, 
and 10 fiducial bundles for acquisition and guiding. Fibres can 
be positioned over a field of radius 30', though there is some vi- 
gnetting starting at r = 20', and image quality outside r — 25' 

1 Throughout this paper, 2MASS identifications are given with coor- 
dinates only, to save space. 

2 http : // cdsweb . u- strasbg . f r/ deni s . html 

3 |http : // www . star . bris.ac . uk/~mbt/topc at/ 1 
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Fig. 2. A sample of G-K low-luminosity supergiants observed 
with AF2. Note the evolution of the Tii 8683A/Fei 8679 A ra- 
tio with spectral type. All these objects have relatively low ex- 
tinction. In most of them (but not in J18392955-0544222) the 
8620A DIB, blended with Fe i 8621 A, is hardly noticeable. 



is quite poor. Due to the size of the fibre buttons in AF2, fibres 
cannot be placed closer than approximately 25" from each other. 
Because of this, it is not an appropriate instrument to observe 
dense clusters. 

We used the R1200R grating, blazed at 7200A. This grating 
provides a nominal dispersion of 0.32 A/pixel over 1300A (this 
is the range common to spectra from all fibres; spectra from dif- 
ferent fibres can be shifted up to about 5% of the range, depend- 
ing on their position in the field). When used in combination 
with the 1'.'6 fibres, it gives a resolution element of 2A. We used 
central wavelength 8300A. Our common spectral range is thus 
7650-8950A, at resolving power R > 4 000. 

AF2 uses as detector a chip mosaic, consisting of two EEV- 
42-80 thinned and AR coated CCDs butted along their long axis 
to provide a 4096 x 4096 pixel detector. As the cameras are 
thinned, they suffer from heavy fringing in the far red spectral 
region (~40%). This effect, however, can be corrected with the 
observation of internal flat-fields. As AF2 is a bench-mounted 
spectrograph, the fringing pattern is very stable. We obtained 
flats at each pointing position to improve fringing removal. Each 
observation was divided into 5-6 exposures. Flat fields were 
taken immediately before and after each of the exposures. With 
this strategy, we have been able to reduce the fringing level 
to %2% for all targets with reasonably bright spectra (^ 4000 
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Fig. 3. A sample of late M-type stars, illustrating the character- 
istics used for classification. The depth of the TiO bandhead at 
8860A is expected to depend linearly with T e e. The 8620A DIB 
is blended with the VO bandhead at 8624A, only visible in stars 
later than M6. The spectrum of the M7.5 1 RSG MY Cep, in the 
open cluster NGC 7419, is shown for comparison. All the spec- 
tra have been normalised in the same, rather arbitrary, way. 



counts/pixel). Except in a few spectra with low count levels, the 
fringing pattern is not visible in the spectra. 

The sky spectra in the far-red are dominated by strong air- 
glow emission lines. For a better sky removal, we used the 
"beam-switch" strategy (moving the telescope back and forth 
between objects and neighbouring sky). For each exposure, we 
moved the telescope ~3' in a random direction and took a sky 
exposure with the same integration time as each of the individ- 
ual target exposures. This procedure was used during the nights 
observed in visitor mode, but not for the service observation. 
Though it is time-consuming, it has allowed a very good sky 
subtraction in most cases. Again, the worst results are found for 
the faintest targets. The observations from June 5th have had the 
sky subtraction carried out with spectra obtained through a high 
number of fibres allocated to blank sky positions. Though the 
procedure is in most cases satisfactory, it gives poorer results 
than the "beam-switch" strategy for the faintest sources. This is 
due to the very different throughputs of the fibres. A few of the 
spectra (with S= 3000 counts/pixel - meaning $ 500 counts/pixel 
in an individual exposure) were unusable due to a combination 
of sky and fringing contamination. 

To observe Stephenson 2, the field of AF2 was centred on 
RA: 18h 40m 25.0s, Dec: -05°44'30". This is not the nomi- 



nal centre of the cluster, but lies ~15' to the Northeast, leav- 
ing the cluster in one of the AF2 quadrant s, and extend i ng the 
field in the direction towards the area where Garzo n et alj (1997) 
had observed an over-density of RSGs. Three configurations 
were observed centred on these coordinates. The first configu- 
ration contained only bright targets, with I < 11.1. The second 
configuration contained targets of intermediate brightness, with 
10.9 < / < 13.1. The third configuration contained faint targets 
with 12.9 < I < 15.1. A small overlap in magnitudes was al- 
lowed, so that a few stars could be picked more than once by 
the configuration software, even though no attempt was made to 
ensure this overlap. 

As a further check, a single configuration with intermedi- 
ate exposure times was observed, displaced a few arcmin further 
Northeast from the main position. The area covered by this con- 
figuration and the three exposures of the Stephenson 2 field rep- 
resent the "main target area" (see Fig. [7]). Finally, to assess the 
singularity of the region surrounding Ste 2, we observed a con- 
trol field, located ~1°1 away from the centre of our main field, 
and so not overlapping with it. The control field was chosen to 
be in an area where no over-density of RSGs had been claimed 
before. In addition, it contains a star cloud, an area with a very 
high density of objects, likely caused by low foreground extinc- 
tion. The densest part of this star c loud is identifiable as the c lus- 
ter candidate Andrews-Lindsay 5 (Andrews & Lindsa vTl967l) . In 
this control region, the selection criteria were loosened in order 
to probe a larger region of the parameter space. We observed a 
few stars with 7.0 < K$ < 7.5, and we allowed a few objects as 
faint as / = 16.1 to be included. 

Due to the small overlap in magnitude between the three ex- 
posures, and the observation of an overlapping field, ~20% of 
the targets have more than one spectrum (with 10 objects having 
three spectra). This has allowed a good estimation of the preci- 
sion of spectral type classifications, of the impact of signal-to- 
noise ratio on the reduction procedure and of the radial velocity 
determination (see Sect. I3.31 l. In total, we have 233 unique stars 
with usable spectra. Of them, 73 lie in the control region, with 
the rest belonging to the main target area. 

In addition, we observed a few sources using the red arm of 
the ISIS double-beam spectrograph, also mounted on the WHT 
The observations were carried out on 2010 June 1st using the 
R1200R grating centred on 8700A. The instrument was fitted 
with the Red+ CCD, which has very low fringing in this range, 
and allows reduction of the fringing to < 1% simply with the 
use of internal lamp flats. This configuration covered the 8380- 
9 100 A range in the unvignetted section of the CCD with a nom- 
inal dispersion of 0.26A/pixel. We used a Of.' 9 slit, resulting in a 
resolving power, measured on the arc lamps, R as 16 000. 

The original purposes of this observation were improving 
the radial velocity determination, assessing the impact of resolu- 
tion on spectral classification and performing a general "sanity 
check" on our results. Unfortunately, the night was far from op- 
timal, with heavy dust cover leading to important extinction and 
poor seeing (> l'/5 most of the time). In consequence, we could 
only observe 1 1 of the AF2 targets, and a few very bright can- 
didates that had not been observed with AF2. Figure Q] displays 
the spectrum of one RSG observed with both instruments so that 
the effects of very different spectral resolutions can be assessed. 

The AF2 data were reduced following the recipe imple- 
mented in the pipeline provided by the AF2 webpag^3 at that 
time. The code follows the standard steps for multi-fibre spec- 



http://www.ing.iac.es/astronomy/instruments/af2/reduction.html 
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trum reduction, using IRAI0 tasks for bias and fiat-field correc- 
tion, throughput calibration between fibres, wavelength calibra- 
tion and sky subtraction. The ISIS data were reduced using stan- 
dard IRAF tasks for bias and flat-field correction, wavelength 
calibration and sky subtraction. 

3. Spectral classification 

The classification of red luminous s t ars in t he far red region 
was discussed by iKeenan & Hvnekl d 1945b . ISharplessI {l956) 



and ISolj (l978) for photographic spectrograms. More recently, 
Paper Jexplored se veral criteria presented in works discussing 
digita l spectra (e.g..lqnestetetalj|i994t ICarquillat etal]|1997t 



iDfaz et al.lll989blZhoulll991h . Here we use the most relevant re- 
sults found to classify our targets. Representative spectra are pre- 
sented in Fig. El EIU and 

3.1. Temperature determination and spectral types 

The far-red spectra of M-type stars are dominated by molecu- 
lar TiO absorption bands, with VO ban ds also becoming promi - 
nent for spectral types M6 and later (Kirkratri ck et al.lll99ll) . 
These bands are much stronger in luminous stars (giant s and su- 
pergiants) than in stars close to the main sequence. In iPaper 3. 
we calibrated the dependence of the depth of the TiO 5(0,0) R- 
branch bandhead at 8860A with T e s, finding a linear relation 



r eff = (3848 + 38) - (1803 + 154) x TiO [K] 



(1) 



valid for giants and supergiants. In this expression, TiO is an 
index defined by 



TiO = 1 - 



^TiO 



(2) 



where the intensities are measured as described in IPaper 1 
As the spectral sequence is a temperature sequence, the intensity 
of this jump is the primary spectral type indicator for M-type 
stars in the range observed. The intensity of the triple TiO band- 
head at 8432^42-52A is an excellent secondary indicator. Some 
examples can be seen in Figures [3] and [4] The TiO bandhead at 
8432A is blended with the atomic doublet Ti i 8435A. For spec- 
tral types earlier than MO, Ti i 8435A is comparable in strength 
to Ti i 8427 A (see Fig. Ell- Between Ml and M3, the former fea- 
ture, blended with the bandhead, grows in strength (see Fig. @}. 
By M4, the ban dhead features at 8442 and 8452A become visi- 
ble dSoiil ll978). By M6, most of the atomic features shortwards 
of 8600A have been obliterated by the growing strength of TiO 
bands, and other bandheads, such as 8558, 8571A are clearly 
seen. By M7, the VO bandheads start to a ppear and the Can 
lines become very weak in giants (Solfl ll978h . as illustrated in 
Fig.S 

The TiO bandhead at 8860A becomes first noticeable, as 
a weak inflection in the continuum at Ml. For earlier-type 
stars, other criteria must be used. The most effective criterion 
is the ratio between the Tii 8683A and the Fei 8679A lines, 
which is very sensitiv e to temperature in the G5-K5 range 
(ICarquillat et al.lll997l) . These two lines are not fully resolved 
at the resolution of the AF2 spectra, but their relative intensities 



5 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



can be discern ed and used to provid e spectral types. Using a ll the 
supergiants inlCenarro et all (12001 1) . ICar quillat et al l (1 19971) and 
iMunari & Tomasellal d!999l) . we find that the ratio is not strongly 
affected by resolution. The Ti i line is very weak in early-G stars, 
and grows after G3 until the two lines have similar depths at G8. 
The Fe i line becomes very weak compared to the Ti i one at K5 
(see Fig. E). This leaves a small range K5-M0.5, where differ- 
ences in the spectral range considered are very subtle. Outside 
this narrow range, spectra can be classified with an accuracy of 
±1 subtypes. 

3.2. Luminosity calibration 

The luminosity of late-type stars can be calibrated with several 
indicators. In general, the intensity of the whole metallic spec- 
trum is dependent on luminosity. The best calibrated indicators 
are the equivalent widths (EW) of the two strongest lines in the 
Can t riplet and the Tii dominated blend at 8468 A ( Diaz et al.l 
119891 iGinestet et al.|[T994llPaper J) , which are useful in the CO- 
M4 and G5-M3 ran ge, respectively. Other lin es sensitive to lu- 
minosity are listed bv lKeenan & Hvnekl d 19451) . They include the 
Fei 8327A line and the blend (mainly Fei) at 8514A. The be- 
haviour of these lines was calibrated in early-M stars, and their 
behaviour at later types has not been calibrated. In any case, due 
to the growing strength of TiO bandheads, they cannot be used 
for types later than M4. It is so necessary to derive a spectral 
type first and then compare to stars of the same type to obtain a 
luminosity class. 

For spectral types later than ~M4, the measurement of most 
features becomes very difficult, as the spectra become domi- 
nated by TiO bands. As a consequence, the luminosity classi- 
fication must be achieved by considering the relative strength 
of the whole metallic spectrum between ~8660A and 8860A. 
This includes, apart from Can 8862A, the complex formed by 
Fei 8675 A, Tii 8683A and Fei 8689 A the 8735A blend, th e 
Mgi triplet at 8807 A and Fei 8824A dAller & Keenanlll95ll) . 
The luminosity class of late M stars with noisy spectra (where 
the strength of weak metallic lines is difficult to asses) cannot be 
considered very accurate. Note also th at there are no MK stan- 
dards for supergiants later than M5 (cf. lSoHf l978). 

In view of the difficulties in classifying late-M stars, we 
have resorted to some approximate classifications. Some stars 
have strong TiO bandheads, indicative of spectral type M7 (i.e., 
TiO ~ 0.5), but show relatively strong metallic lines and no sign 
of the VO 8624A bandhead. We have classified these objects as 
M6-7. Many other stars show evidence for the VO bandhead, 
and are therefore M7 or later, but their spectra are rather noisy 
(likely because the airglow is more difficult to remove due to 
much lower count rates inside the huge TiO features). In some 
cases, the metallic lines are quite strong, clearly indicating that 
the star is not a giant, but we cannot be sure that they are suf- 
ficiently strong to be a supergiant. In such cases, we have used 
the I— II classification, which is meant to indicate that the star is 
likely a very late supergiant, but may also be a bright giant. 

3.3. Radial velocities 

Radial ve locities were calculated using the same algorithm de- 
scribed in Paper I. In that paper, we cal culated cross-correla tions 
against the spectra in the library of ICenarro et al.l d200ll) . In 
this work, however, we cross-correlate our spectr a against a set 
of m odel spectra from the POLLUX database (Palacios et al. 
2010). This change has been motivated by two reasons: firstly, 
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Fig. 4. A sample of low-luminosity M-type supergiants observed 
with AF2. All these objects are likely foreground to Ste 2, as 
they have lower vlsr and lower E(J - K§) colour excess. In 
these sources, Fei 8621 A is stronger than Fei 8611 A, indicat- 
ing an important contamination by the 8620A DIB. The main 
luminosity indicators are the Ca n lines and the 8648 A blend, 
but note also the weakness of the 8514A feature compared to 
objects in Fig. [5] The growing strength of the TiO 8432A band- 
head, blended with Ti 1 8435 A, can be seen to correlate very well 
with TiO 8660A (the Ti 1 8427 A line, immediately to its left, can 
be used as reference). 



the higher resolution available for the synthetic spectra makes 
them more suitable for the dataset we use here (especially, for 
the spectra obtained with ISIS , with R > 15 00). S econdly, the 
radial velocity corrections of ICenarro et al. I d200l are not al- 
ways accurate. Many luminous red stars - particularly those of 
later spectral types - displa y radial veloci t y varia tions, but most 
of the stars in the library of ICenarro et al.l Q2001 ) were observed 
only once and then corrected with a shift reflecting their aver- 
age radial velocity. As a consequence, the use of these spectra as 
templates introduces additional dispersion in radial velocities. 
By correlating all the spectra in the library against the POLLUX 
synthetic data, we have found that this dispersion can be as high 
as 8 km s _1 for some given stars. 

The synthetic spectra were degraded to the instrumental res- 
olution of our spectra. To reduce the uncertainty in radial veloc- 
ity due to the choice of a template, every spectrum was calibrated 
against its closes t match in spectral type, following the proce- 
dure described in Paper I. The internal accuracy of our measure- 
ments can be gauged from the several stars that were observed 



Fig. 5. A sample of high-luminosity supergiants observed with 
AF2. All these objects suffer moderately high to high extinc- 
tion, as indicated by the 8620A DIB blended with Fei 8621 A. 
All these objects have radial velocities similar to members of 
Ste 2, except 18395861-0549549, which has a slightly lower 
Vmd- The strength of the 8514A feature compared to the neigh- 
bouring Ti i 85 1 8 A is an indicator of high luminosity. 



two or three times with AF2. The differences between the spec- 
tra were measured and the dispersion in v ra d was found to be 
<x(v) = 2.8 kms -1 . This value is surprisingly good when com- 
pared to the spectral resolution, but the fact that AF2 is a bench 
mounted spectrograph results in very high stability. 

Our two sets of spectra (from AF2 and ISIS) were observed 
along with stars of known radial velocity, allowing us to derive 
absolute errors for our measur emen ts. The AF2 data include 3 
supergiants already observed in D07. The average difference be- 
tween our values and those of lD07l is Av rat j = — 1 ± 4 kms~', 
indicating that there are no systematic differences. 

For the ISIS dataset, 5 red giants wit h highly accurate v ra d 
measurements in iMermilliod et al.l d2008l) were observed along 
w ith the science targets. Comparing our measured v rac j to those 
of Mermil liod et al.l ( 20081) . we find a difference Av ra( j = — 1.2 ± 
0.3 kms -1 . Such difference is compatible with the highest accu- 
racy in radial velocity determination achievable with ISIS, ac- 
cording to its manual. 

Finally, we can compare the radial velocities in our two 
datasets. The ISIS and AF2 datasets have 1 1 stars in common. 
Once the systematics have been removed, their velocities agree 
very well, with Av rac j = 1+5 kms~'. In view of all this, we 
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estimate that our radial velocities have a conservative total un- 
certainty of +5 km s -1 . 

The radial velocities measured were transformed into the 
Local Standard of Rest (LSR) reference system using IRAF's 
rvcorrect package. In this direction (7 » 24° - 30°), the Galactic 
rotation curve allows the separation of stars located in differ- 
ent Galactic arms . The distribution of molecular clouds along 
the line of sight dRathborne et al.l 120091) shows a small group 
of objects with v'lsr~ +25 kms" 1 , perhaps identifiable as the 
Sagittarius Arm (but see below). Then we find a number of 
sources with velocities in the -40-65 kms -1 range, generally 
interpreted as a first crossing of the Scutum-Crux arm. There is a 
clear gap for v LSR 70 - 90 km s -1 (except for the 28° - 30° 
range), after which we find a high number of sources with ve- 
locities between 90 and 115 km s -1 . lD07l determine an average 
systemic vel ocity vlsr = 1 10k m s _1 for Ste 2. The distribution of 
H ii regions dBania et al.ll2010 ) confirms the low number of trac- 
ers associated with the Sagittarius arm, but reflects a more ho- 
mogeneous distribution in velocities between 40 and 90 kms -1 , 
with again a high concentration in the vlsr ~ 90-115 kms -1 
range. 

4. Results 

Once again, the selection criteria have been very effective. All 
the objects observed are late-type luminous stars, with one single 
exception, J18390349-0546386, which is a B-type supergiant0. 
The rest of the AF2 targets comprise 3 G-type supergiants, a 
number of K-type stars with luminosity classes II and higher 
and M-type giants and supergiants. Early M-type stars have lu- 
minosity classes II and higher, while for spectral types later than 
M3, we also detect luminosity class III objects. This distribution 
of spe ctral types seems to confirm the hypothesis advanced in 
Paper I that our 2ir cut selects stars with very extended atmo- 
spheres. Table Q] presents the derived spectral types for all our 
targets and th e colour excesses E(J - K$), calculated using the 
calibration of iLevesque et al.l d2005h . 

Table |2]presents the same data for the ISIS targets. There are 
eleven sources in common with the AF2 sample, and spectral 
types and vlsr have been derived independently with the aim 
of gauging the effect of resolution on the procedures used. The 
differences in spectral type are minimal, with an occasional very 
small disagreement in the luminosity class. 

The number of true supergiants (i.e, objects with spectral 
types indicative of high-mass starsQ) detected is quite high. The 
AF2 sample includes more than 60 supergiants (72 if we assume 
that late-M stars classified I-II and with vlsr > 90 kms -1 must 
be supergiants) and 1 1 others with Ib-II classification. The ISIS 
sample adds 6 more RSGs. Of these 78 RSGs, only five are mem- 
bers of Ste 2 known from previous works. As a comparison, be- 
fore the discovery of the RSG clusters, the number of confirmed 
RSGs in the Galaxy was just a few hundred. 

We find an unexpectedly high number of stars with veloc- 
ities in excess of vlsr = +120 kms -1 , i.e., higher than any 
value expected along this line of sight, with a few values vlsr ^ 
+ 145 kms -1 . Such stars must have high peculiar velocities. 

The presence of stars with very high peculiar velocities 
is demonstrated by several measurements of negative vlsr, 
which are not expected at all along this line of sight. Most 

6 The 2MASS magnitudes of this object clearly correspond to a late- 
type star. The DENIS J and K magnitudes are compatible within the 
errors. The object is saturated in UKIDSS, but, with the higher spatial 
resolution, it is clearly resolved into two stars, separated by ~ 1'.'8. We 
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Fig. 6. Spectra of five RSGs in Stephenson 2 and one possible 
runaway member (D39). The top four spectra are from AF2, 
while the bottom two are from ISIS, degraded to a similar res- 
olution. Star D2 is a very late RSG, associated with an IRAS 
source and maser sources (compare to MY Cep in Fig. [3). The 
spectrum of D39 has very low quality, as this is a faint target. 



of these objects are late M giants, the most extreme case 
being J 18403943 -05 16043 (M4.5III) which displays v LSR = 
-72 kms" 1 . This object presents some excess in the far-IR (the 
21 -/mi WISE point deviates considerably from a stellar model), 
but the 5/?/fzer/GLIMPSE data do not indicate mid-IR excess. 
Therefore most of the E(J - K&) w 0.65 must be interstel- 
lar and this should be a relatively distant object, located no 
closer than the Sagittarius ar m. Its peculiar velocity component 
is thus at least -100 kms -1 . Deg uchi et al.l (1201 Obi) found sev- 
eral maser sources with similar velocities at Galactic longitudes 
I = 20° - 40°, and interpreted them as members of streaming 
motion groups. Stars with vlsr ~ 120 km s , on the other hand, 
might be following orbits profoundly disturbed by the Galactic 
Bar potential. Most of them have moderately high colour excess, 
suggesting that they are distant objects. As these objects might 
be signposts of systematic deviations from the Galactic rotation 
curve, the accuracy of kinematic distances along this line of sight 
may well be rather lower than generally assumed. 



believe that the /-band magnitude is likely dominated by the second 
star, as there is no evidence for a late-type component in our spectrum. 

7 These would be G and K stars of luminosity classes la and lab, and 
all M-type supergiants. 
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Table 1. List of AF2 targets with usable spectra with their main properties. 



ID 




./ 


H 




E(J - K s ) 


Vlsr 


Spectral 














(kms ') 


type 


18382013- 


-0543434" 


8.37+0.02 


6.70 + 0.04 


5.95 ±0.01 


1.35 


+59 


Ml. 5 lab 


18383438- 


-0527029 


8.94 ± 0.05 


7.52 ± 0.02 


6.95 ± 0.01 


0.70 


+67 


M7 III 


18384011- 


-0547537 


7.42 ± 0.01 


6.10 ±0.03 


5.48 ±0.01 


0.64 


+59 


M7 II 


18384958- 


-0557222 


8.95 ± 0.01 


7.20 ± 0.03 


6.39 ± 0.02 


1.27 


-12 


M7 III 


18385040- 


-0549183 


9.81 +0.02 


7.78 ± 0.04 


6.86 ± 0.03 


1.70 


+91 


M6III 


18385144- 


-0600228 (D12) 


7.22 ± 0.04 


5.94 ± 0.02 


5.35 ±0.01 


0.62 


+35 


M7 II-III 


18385268- 


-0544207 


8.30 ± 0.01 


7.08 ± 0.02 


6.62 ± 0.01 


0.62 


+43 


M1.5Ib-II 


18385336- 


-0557444 (D38) 


8.57 ±0.01 


7.32 ± 0.01 


6.79 ± 0.03 


0.61 


-4 


M4.5 II 



8 The JHKs magnitudes are from 2MASS. Spectral types are b ased solely on our sp ectra. The colour excess has been calculated from the 
observed (J - K s ) colour and the tabulated intrinsic colour, after lLevesque et al] (l2005). 

Table 2. List of ISIS targets. Objects marked with a t were not primary targets and their identifications are based on slit reconstruc- 
tion. Such objects felt in the slit by chance and do not necessarily fulfill the selection criteria. 



ID 




./ 


H 




E( J - K s ) 


Vlsr 
(km s" 1 ) 


Spectral 
type 


18365544- 


-0643364" 


8.35 ± 0.02 


7.34 


±0.02 


6.95 ± 0.03 


0.51 


+48 


K4II 


18372648- 


-0628397+ 


8.88 ±0.03 


6.95 


±0.02 


6.09 ± 0.02 


1.73 


+68 


Ml lb 


18372699- 


-0628444 


7.12 ±0.02 


5.90 


±0.02 


5.37 ± 0.02 


0.50 


+66 


M6-7 II 


18382442- 


-0522521 


8.14 ±0.03 


7.00 


±0.03 


6.57 ± 0.02 


0.44 


-21 


M3.5 II-III 


18383565- 


-0603043 


7.78 ± 0.02 


6.64 


±0.02 


6.18 ±0.02 


0.47 


+36 


M3.5 III 


18384011- 


-0547537 


7.42 ± 0.02 


6.10 


±0.03 


5.48 ± 0.02 


0.64 


+24 


M7 II-III 


18391838- 


-0600383 (D6) 


7.72 ± 0.02 


5.92 


±0.03 


5.07 ± 0.02 


1.51 


+ 108 


M3.5Iab 


18391952- 


-0559194 1 " (D19) 


8.28 ± 0.02 


6.58 


±0.04 


5.80 ± 0.02 


1.38 


+ 111 


M2 la 


18391961- 


-0600408 (D2) 


6.90 ± 0.02 


5.05 


±0.02 


4.12 ±0.27 


1.49 


+ 121 


M7.5I 


18392506- 


-0610036 


7.76 ± 0.05 


6.53 


±0.04 


6.05 ± 0.02 


0.58 


+61 


M3 II-III 


18394635- 


-0559473 


7.05 ± 0.02 


5.34 


±0.02 


4.51 ±0.02 


1.36 


+ 110 


M4Ib 


18395282- 


-0535172 


6.81 ±0.02 


5.41 


±0.03 


4.75 ± 0.03 


0.74 


+56 


M8 1-II 


18402741- 


-0556336+ 


12.09 ± 0.03 


10.84 


±0.02 


10.38 ±0.03 


0.58 


+21 


M3 II 


18402744- 


-0556176 


8.58 ± 0.02 


7.44 


±0.02 


6.95 ± 0.02 


0.50 


+53 


M3 II 


18403179- 


-0529522 


7.88 ±0.03 


6.55 


±0.04 


5.98 ± 0.02 


0.61 


+36 


M7 II-III 


18405021- 


-0552520 


7.66 ± 0.03 


6.49 


±0.03 


6.00 ± 0.02 


0.41 


+71 


M6 II-III 


18405351- 


-0517488 


8.42 ± 0.03 


7.20 


±0.04 


6.74 ± 0.03 


0.69 


+71 


K7Iab 


18405360- 


-0525004 


7.53 ±0.02 


6.12 


±0.05 


5.46 ± 0.02 


0.93 


+ 113 


M3Ib 


18410764- 


-0535045 


8.16 ±0.02 


7.12 


±0.03 


6.67 ± 0.02 


0.25 


-26 


M6-7 II 


18410851- 


-0555046 


8.18 ±0.02 


6.98 


±0.03 


6.50 ± 0.02 


0.51 


+71 


M4.5 II 


18412383- 


-0526073 


7.79 ± 0.02 


6.43 


±0.03 


5.87 ± 0.02 


0.90 


+ 103 


M0 lab 


18414143- 


-0531378 


8.23 ± 0.04 


7.23 


±0.04 


6.81 ±0.02 


0.61 


+73 


K3 lab 


184151 12- 


-0610503 


7.99 ± 0.02 


6.91 


±0.03 


6.48 ± 0.03 


0.31 


+21 


M5 Ib-II 


18415256- 


-0611381+ 


11.83 ±0.02 


10.13 


±0.02 


9.42 ± 0.02 


1.78 


-1 


K0III 


18421965- 


-0545389+ 


10.32 ± 0.02 


9.30 


±0.03 


8.87 ± 0.03 


0.44 


+49 


M0 lab 


18422147- 


-0545383 


8.53 ± 0.02 


7.41 


±0.05 


6.98 ± 0.02 


0.37 


+34 


M4.5 III 



8 The JHK S magnitudes are from 2MASS. Spectral types are b ased solely on our sp ectra. The colour excess has been calculated from the 
observed (J - Ks) colour and the tabulated intrinsic colour, after Leves que et alj J2005). 



Table 3. New spectral types for members of Stephenson 2. 



ID 


Co-ordinates (J2000) 




./ 


H 


Ks 


Spectral 


Spectral 


Vlsr 


vlsr D07 




RA 


Decln. 










Type (CO) 6 


Type 


(kms" 1 ) 


(kms" 1 ) 


D2 


18 39 19.6 


-06 00 40.8 


12.66 


6.90 


5.05 


4.18 


M3I 


M7.5I 


+ 121 


+ 111 


D6 


18 39 18.4 


-06 00 38.4 


12.65 


7.72 


5.92 


5.06 


M5I 


M3.5 lab 


+ 108 


+ 107 


D15 


18 39 22.4 


-06 01 50.1 


13.07 


8.13 


6.35 


5.51 


M2I 


M1.5Iab 


+ 106 


+ 112 


D18 


18 39 22.5 


-06 00 08.4 


12.88 


8.18 


6.45 


5.63 


M4I 


M0.5 lab 


+ 114 


+ 111 


D19 


18 39 19.5 


-05 59 19.4 


13.01 


8.28 


6.58 


5.80 


M2I 


Ml lab 


+111 


+ 107 



8 JHKs magnitudes are from 2MASS. i magnitudes from DENIS. 
b Estimated from the strength of the CO 2.3 pm bandhead bv lD07l 
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4.1. Selection biases 

To assess the statistical value of our results, we need to evalu- 
ate the importance of selection biases in the sample finally ob- 
served. The most obvious source of bias is the fact that we are 
observing highly-reddened sources (almost) in the optical range. 
The choice of the observation range is determined by the high 
multiplexing capability of optical spectrographs, but also by the 
much better ability to obtain spectral classifications for red lumi- 
nous stars in the 8300-8900A range than in the K band, where 
our targets are intrinsically much brighter and extinction is much 
lower. 

Because of this, it is encouraging to find that a very sub- 
stantial fraction of the objects passing the 2MASS cuts, have a 
counterpart in the / band (note that the i or / bands fall on the 
short wavelength side of our spectral range). In the main target 
area, 73% of the 2MASS candidates have an /-band counterpart 
with / < 15.1, and where so eligible for observing. In the con- 
trol zone, the fraction was 68% (though a few targets with fainter 
counterparts were also observed). Therefore the fraction of can- 
didates too faint to be observed with our instrumentation was 
~ 30%. 

Finally, it is worthwhile mentioning that there is no clear cor- 
relation between brightness in the / band and spectral type. Even 
though a significant fraction of the catalogued members of Ste 2 
have / magnitudes that would have put them close of the faint 
end of our "intermediate-brightness" configuration, we find a 
very significant number of supergiants in the "bright" configu- 
ration. Many of them are clearly foreground to Ste 2, but there 
are a few RSGs with vlsr > 90 km s~' and rather high reddening 
that still manage to have / < 1 1 . 1 . In view of this, our "bright" 
cut at / < 10 may perhaps seem too restrictive. In any case, less 
than 1% of the 2MASS candidates were rejected by this magni- 
tude cut. 

The Ks = 7.0 cut was chosen as ~90% of the confirmed 
members of Ste 2 are brighter than Ks = 6.2. If a given star is 
fainter than the cluster members because of higher extinction, 
one would expect its / magnitude to be much fainter than those 
of cluster members, making it an unpractical target for this sur- 
vey. If, on the other hand, a star is fainter because of intrinsic 
lower luminosity, it is likely to represent an older or more dis- 
tant population. Setting the limit at Ks =7.0 seemed a sensible 
approach to account for these possibilities, while still keeping a 
reasonably small sample. 

Simply moving the magnitude cut to Ks = 7.5 increases 
the total sample (objects observable in / plus those not acces- 
sible) from slightly below 600 to more than 900 stars. For this 
increased sample, the fraction of candidates with a detection in 
/ and the fraction of candidates with 10 < / < 15.1 are not 
significantly different from the Ks = 7.0 sample. Analysis of 
the 2MASS CMDs and Q/K s diagram, however, shows that for 
A"s ^ 8 field contamination becomes so important that a sub- 
stantial number of objects fulfilling our colour criteria could be 
foreground F-type stars. 

4.2. Stephenson 2 

Care was taken to include some members of Ste 2, with ra- 
dial velocities measured by lD07l to be used as calibration. 
Unfortunately, it is very difficult to place AF2 fibres on mem- 
bers of compact clusters, because they collide. Three stars were 
observed with AF2, one of them (D19) included in two config- 
urations. D19 and two additional stars were observed with ISIS. 
The list of Ste 2 stars is given in Table [3j together with the de- 



rived spectral types and a comparison of their properties in lD07l 
including the spectral type estimation based on the strength of 
the CO 2.3 /im bandhead. The spectra are displayed in Fig. [6] 

We observed two stars in common with Paper I. The spectral 
type of D15 is Ml. 5 lab, as in the previous paper. The spectral 
type of D18, M0. 5 lab, is slightly different from that found in 
Paper I (Ml. 5 lab), though we must note that differences in this 
spectral range are v ery subtle. This star was classified as M4 1 by 
Davi es et~al] ([2007). Another star presenting a large difference 
in spectral type is D2, whichlD07l give as M3 1, but appears very 
late here (we give M7.5I). This new classification, when con- 
sidered together with its identification with the mid-IR sources 
IRAS 18366-0603 and MSX6C G026.2038-00.0574, strongly 
points to this star (the brightest confi rmed cluster member in Ks) 
being surrounded by a dust envelope. Degu chi et alJ (1201 Oal) find 
both SiO and H2O masers associated with this object (see later). 

Our sample includes a few candidates from [DOl that were 
considered non-members or did not have their radial veloci- 
ties measured. The first one is J18385144-0600228 (D12). Its 
moderate luminosity (M7II-III) and its vlsr = +35 kms~' rule 
out membership. Its low reddening E(J - A" s ) * 0.60 is typi- 
cal of the first crossing of the Scutum-Crux arm, at ~ 3.5 kpc. 
Likewise, J 18385336-0557444 (D38) is a bright giant (M4.5 II) 
with vlsr = -4km s _1 . Again, its low reddening E(J-K S ) ~ 0.60 
confirms that it is a foreground star. 

J 1839 1636-0603 149 (D28) was considered a foreground gi- 
ant by lD07l because of the low EW of the A"-band CO bandhead. 
Our spectrum (see Fig. [2}, however, shows that it is a super- 
giant. The CO bandhead is weak because of its early spectral 
type, G5 lab. This star is projected very close to the central clus- 
ter concentration. Its vlsr = +27 kms~' rules out, in principle, 
a connection to the cluster. Its E(J - Ks) = 1.0 is lower than 
that of cluster members and the 8620A DIB is rather weak in its 
spectrum. This object is thus foreground to the cluster. 

Finally, J18390701-0605043 (D39) is an Ml lab supergiant. 

Its 

vlsr — +63 km s is very different from those of clus- 
ter members, but its E(J - Ks) = 1.79 is typical for cluster 
members^. This star lies between Dl a nd D5, in the Southern 
concentration that lDeguchi et alJ (l2010al) call Ste2 SW. Its spec- 
tral type is also typical of cluster members. If the star is fore- 
ground, it must be rather less luminous than its luminosity class 
indicates. Because of this, we cannot rule out the possibility that 
this is a cluster member with peculiar velocity. 

Co nverse ly, we find several stars outside the region stud- 
ied by |D07| that are very likely connected to the cluster. For 
example, J 18393774-0556 183 is an M0 la supergiant with 
E(J - K s ) = 1.39 and v LS r = +98 kms" 1 only ~ 7' to the 
NE of the centre of Ste 2. Two arcmin further to the North, 
we find J18392977-0553032 (StRS 236), an M3.5Iab super- 
giant with E(J - K s ) = 1.48 and v LSR = +114 kms~'. 
Slightly to the South, and just outside the region scanned 
bv lD07l we find J18394635-0559473 (IRAS 18370-0602 = 
MSX6C G026.26 78-00.1492), M4.5 1 already classified as a 
candidate RSG bv lRawlings et al.l(l2000l) . lDeguchi et al.ld2010al) 
excluded an association of this object with Ste 2 due to its 
spatial distance to the cluster, but the measured radial velocity 
vlsr = +114 kms -1 and E(J - Ks) = 1.36 are comparable to 
those of cluster members. 

Continuing in this area, J18393077-0551103 is an Ml lab 
supergiant ~11' North of the cluster core, with very high 



8 The median value for the clust er is E (J - A"s) = 1 .76, calculated for 
the 26 likely members selected by D07, using the same calibrations as 
in that work. 
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reddening E(J - K s ) = 1.98 and v LSR = +114 kms -1 . 
J 18393544-0548 186 is a K7Iab supergiant -14' away from 
Ste 2, with E(J - K s ) = 1.23 and v LSR = +106 kms" 1 . 
Slightly further away, we find J18394562-0544057, M4.5Ib, 
with£'(/-A's) = 1.08 and vlsr = +H3kms _1 . All these objects 
seem to define an extension of the cluster towards the Northeast, 
without a clear divide. 

Interestingly, we also find two late-M giants with similar 
velocities, J18393324-0529435, M7.5III (J - K s = 3.2), and 
J18393145-0537096, M7III (J - K s = 2.8). Both are very 
red objects, and an important part of the obscuration must be 
interstellar^. In view of this, they are likely to be at the same 
distance as Ste 2. We find several other examples of late M gi- 
ants with LSR velocities > 90 kms 1 and high reddenings at 
larger angular separations from Ste 2. 

4.3. Spatial distribution of RSGs 

Figure [7] shows the spatial distribution of the objects classified 
as RSGs, compared to the distribution of other sp iral tracers, H n 
region s observed in radio from the catalogue of lAnderson et al.l 
d201 1 ), and methanol m asers (associated to young stellar ob- 
jects; Green et ai1l201 ll) . The two separate circular regions cor- 
responding to the main target area (East) and the control zone 
(West) are indicated. Ste 2 lies on the Eastern side of the main 
target area. The number of RSGs in the circle including Ste 2 
is obviously higher (about two thirds of the total number), but 
this area was observed with four different configurations (plus 
the extra ISIS targets), against only one configuration for the 
control region. In view of this, we find no compelling evidence 
suggesting that RSGs are more numerous in the region surround- 
ing Ste 2, outside the central concentration, than in the control 
region. 

We only observed 71 stars in the control area, out of ~ 156 
possible targets, detecting 27 RSGs (38%). In the main target 
area, however, we observed 162 stars out of 246 possible tar- 
gets, detecting 45 RSGs (28%). Therefore, even taking into ac- 
count the fact that we cannot probe effectively the central region 
of Ste 2, the fraction of RSGs in the control zone is not signifi- 
cantly different from that in the main target area. Moreover, most 
of the 27 RSGs detected in the control zone have vlsr velocities 
> 90kms~'. This is more easily seen in Fig. [8] which shows the 
distribution in vlsr of all the RSGs detected in our survey com- 
pared to other spiral tracers. Over the whole I = 26° - 28° range, 
there is a high concentration of sources with vlsr ~ 115 km s _1 . 
A few methanol masers have similar velocities, while Hn re- 
gions prefer velocities slightly below 110 kms~'. Many other 
RSGs, have velocities between 1 00 and 1 1 0km s~ 1 . This concen- 
tration of sources delimits the terminal velocity in this direction. 
Different Galactic rotation curves imply terminal velocities be- 
twen vlsr = 105 and 120 kms _I , depending on adopted values 
for t he Sun's orbital velocity and distance to the Galactic centre 
(e.g., IClemensIl 19851; iLevine etal . 2008: lReid etZll2009l) . 



9 Theoretical magnitudes calculated from PHOENIX models and cal- 
ibrated against observational data su ggest that all late- M giants have 
(J-K ) ~ 1.3 ([Kucinskas et al. 2005). The calibration of Groenewegen 
d2004h . based on older models, gives (J - K) = 1.37 for dust- 
free M7 stars. This reference also gives (J - K) = 1.97 for dust- 
enshrouded M7 stars. None of the two stars is associated to an IRAS 
source, suggesting that they are not dust enshrouded. We note that 
IRAS 18368-0532 is unlikely to be associated to J18393324-0529435, 
but rather to J 18393238-0530092, which is ~ 30" away and ~ 1.5 mag 
brighter in all GLIMPSE bands. 



We find many RSGs close to Ste 2 with vlsr almost iden- 
tical to cluster members, but a substantial population of RSGs 
with similar velocities extends to very large angula r distances 
from t he cluster. This result confirms the claim bv lGarzon et al.l 
(1 19971) of a strong over-density of RSGs towards / = 27°, and 
demonstrates th at this over-density exists over a much larger area 
than sampled bv lLopez-Corredoira et al.l (Il999l) . 

On the other hand, the most unexpected result of our survey 
is the large number of RSGs with velocities significantly smaller 
than the terminal velocity. We find ~25 RSGs (not counting ob- 
ject classified Ib-II) with velocities vlsr < 90 kms -1 . A few 
of these objects have spectral types compatible with being (but 
not necessarily implying) intermediate-mass (~6-8M0) stars, 
such as J18395686-0603482 (G8Ib) or J18410523-0525076 
(K2 lb). Most of them, however, have spectral types suggestive 
of true supergiants (i.e., high-mass stars). Their number seems 
surprisingly high when compared to the census of RSGs known 
in the Galaxy or the Magellanic clouds, as they are simply the 
foreground population, likely to represent mostly the field popu- 
lation of the Scutum-Crux arm (see below). 

5. Discussion 

We have found high numbers of RSGs, with very different ra- 
dial velocities and reddenings, spread around the galactic equa- 
tor between longitudes / = 26° and / = 28°. These objects must 
be tracing recent events of star formation. Here we try to un- 
derstand their 3D spatial distribution by making use of the data 
available. 

5.1. The sightline to Stephenson 2 

Figure [9] shows the run of the extinction towards three direc- 
tions covered by our survey , calculated using the technique of 
ICabrera-Lavers et al.l (I2005I) on UKIDSS data for circles of ra- 
dius 15'. This technique identifies the red clump giant population 
and estimates the average colour excess at a given distance mod- 
ulus. The three figures are very similar to the c orresponding dia- 
gram s in the directions to Alicante 8 (/ = 24° 6; iNegue ruela et al. 
2010) and the RSGC3 association (/ * 29°2: lPaper J) . Along all 
the sightlines, there is a relatively abrupt increase in the absorp- 
tion around d « 3.5 kpc, identifiable as the first crossing of the 
Scutum-Crux arm, where E(J - Ks) climbs up from low values 
to ss 0.8. After this, the extinction grows gradually until ~6 kpc, 
where it suddenly becomes so high as to render the technique un- 
usable. The line of sight to Alicante 8 is slightly different, as the 
extinction remains relatively stable after 3.5 kpc and the sudden 
increase happens at ~5 kpc, in front of the cluster. 

The run of extinction to Ste 2 is also slightly different in one 
important aspect: the sudden rise in extinction happens slightly 
beyond 5 kpc. Taken at face value, the extinction curve places 
Ste 2 at « 5.3 kpc. This val ue is compatible within th e errors 
with previous determinations dOrtolani et al.l2002l lD07h . though 
it lies near the low edge and forces the age value towards the 
high-end of the permitted range (~ 20 Myr). This distance value, 
though, should not be considered very accurate, as the extinc- 
tion is determined over a 15' circle and there is clear evidence 
of very patchy absorption over the field. Maps of CO emission 
in thi s direction from the Galactic Ring Survey d Jackson et al.l 
2006) show very little emission below v L sr = +90 km s 1 and 
very strong emission for higher velocities. 

All the directions show a very low amount of extinction 
associated with the Sagittarius arm, at ~ 1.5-2 kpc. Our sam- 
ple contains very few possible tracers of this arm, but this is 
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Fig. 7. Spatial distribution of red supergiants detected in our survey. Left panel: Comparison to other spiral tracers. Black filled 
circles represent RSGs observed in our survey. Their size is proportional to the estimated E(J - as shown in the l egend. Open 
circles are other RSGs from ID07I Diamonds are methan ol masers (associate d with star-forming regions) taken from lGreen et al.l 
(I2011h . Asterisks represent H n regions observed in radio d Anderson et al.OOllI) . Right panel: Velocity distribution of RSGs. The 
size of the symbols is proportional to their vlsr- The blue circle represents an RSG with vlsr < . In both panels, the nominal 
position of Ste 2 is indicated by a circle of radius 6'. The large red ovals indicate the field of view of AF2 for the main target area 
(and the overlapping check exposure) and the control area (to the left). 



fully expectable, as our colour cuts are chosen to leave out stars 
with little obscuration. F or example, the open clu ster NGC 6664 
(d ~ 1.4 kpc, / = 24°0: iMermilliod et al.ll 19871 |2008b contains 
several G-K II stars, all of which have I <9 and would thus have 
been rejected by our /-band brightness cut. 

In any case, spiral tracers for the Sagittarius Arm are scarce 
in this area. The maps of molec ular clouds show a few ob- 
jects with velocities «+25 kms -1 (Rathbo rne et al.l l 2.Q09b - The 
Hn regions G026.103-00.069aB (v L sr = +29.3 kmr 1 ) and 
G026.600-00.020a (v LSR = +21.6 kms -1 ) have velocities com- 
patible with this arm, but the fact that they are not visible in 
Ha casts doubt on the possibility that their distance is < 2 kpc. 
Moreover, no clusters with ages < 20 Myr and distances com- 
patible with the Sagittarius arm are known. 

In our sample, there are a few bright giants with velocities 
in this range and relatively low reddenings. Good candidates 
to be located in the Sagittarius arm are J18395177-0549409 
(M5II-III, v LS R = +19 kms" 1 ), which has E(J - K s ) « 0.4, 
or J18392533-0526088 (M2Ib-II, v LSR = +20kms 1 ), with the 
same E(J - The only supergiants that could be located 
in the Sagittarius arm are J18421727-0548553 (M5Ib, v LSR = 
+28 kms- 1 ), with E(J - K s ) = 0.44, and J18391636-0603149 
(= D28; G5Iab, v LSR = +27 kms _1 )El 



10 Radial velocities for Hn regions are from Ande rson et al.l d201lh 
and references therein. 

11 In principle, the observed E(J - K s ) at 1.0 suggests a much higher 
distance. We can estimate the luminosity using the calibration of the 
strength of the Oi 7774A by Arellano Ferro et al. (2003). We measure 
EW= -0.6A, corresponding to M v = -3.0, according to this calibra- 
tion. Using the intrinsic colour for a G5 supergiant, Mk ~ -4.8. This 
would imply a distance of only ~ 1.1 kpc. Of course, the luminosity 
may be underestimated because the strength of the feature is decreasing 
for late-G spectral types (it has disappeared by K0), while the calibra- 
tion has very large uncertainties (standard errors +0.4 mag). However, 
given these values, the source is likely to be in the Sagittarius arm. With 
a G5 lab spectral type, the star may be a Cepheid variable. If this is the 



Towards / = 28 - 29°. lTurnerl(ll980l) found a population of 
O-type stars and B supergiants at d * 3.5 kpc and relatively low 
reddenings, E(B - V) = 1.1-1.3 (similar to the colour excess 
of stars in the Sagi ttarius arm in the same dire ction). Likewise, 
towards / = 25°3, iMarco & N egueruelal (1201 ll) found OB stars 
with E(B — V) at 1.0, and evidence for obscuring clouds at a 
distance ~ 3 kpc. Several Hn regions with v'lsr~ +40 kms -1 
can be found in both directions (see Fig. [HJ. We can identify this 
populations with the near side of the Scutum-Crux arm, just in 
front of the clouds causing absorption at ~3.5 kpc. 

In view of this, we would expect to see stars with vlsr ~ 
40 kms -1 and E{J - K&) w 0.6 tracing the near edge of the 
Scutum-Crux arm in our region. However, we find very few ob- 
jects with such parameters between / = 26° and / = 28°{3 This 
low number of tracers is not due to our colour cuts, as our lower 
limit (J-Ks) > 1.3 will only leave out stars with intrinsic colour 
(J - Ks)o < 0.7, i.e, earlier than ~K0, at this reddening. 

In contrast, between I » 26° and 27°, we find RSGs with 
vlsr = 60-70kms~' and moderately high colour excesses. T hey 
seem to be accompanied by a large number of H n regionsPT 
This is clearly the only structure visible in Fig. [8]in this I range. 

case, its 2MASS colour might not correspond to the observed spectral 
type. 

12 On the Western side, only J18385268-0544207 (Ml. 5 Ib-II, v LSR = 
+43 kms" 1 ), J18390034-0529144 (K4 Ib-II, v LSR = +43 kms-') and 
perhaps J 18385 144-0600228 (M6I-II, v LSR = +35 kms" 1 ) seem possi- 
ble tracers of such structure. These three objects lie quite close together 
and very near the edge of the region surveyed (i.e., close to / = 26°). 
Towards / = 28° we also find a few stars with vlsr ~ +36 kms -1 , but 
much higher reddenings. Objects like 118423886-0446094, Ml. 5 lab 
and E(J - K s ) ~ 1.3, or 118442147-0434460, K2Iab and E(J - K s ) ~ 
1.0, might be RSGs associated with high-obscuration regions in the 
Scutum-Crux arm. 

13 In the catalogue of And erson et al.l d201lh . we find 
G025.469-0.121b (v LSR =60.1 kms" 1 ), G025.477+0.040a 
(v LSR =60.5 kms" 1 ), G26.317-0.012b (v LSR =62.9 kms 4 ), 
G026.600-00.020c (v LSR =69 kms- 1 ), G026.600-00.106 
(v LSR =69 kms- 1 ), G26.610-0.212b (v LSR =69.9 kms" 1 ), 
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Fig. 8. Diagram showing the (vlsr, distribution for different spiral tracers in the vicinity of Stephenson 2. Black circles represent 
red supergiants observed in our survey. Their size is propo rtional to the estim ated E( J - Ks), as in Fig. [7] Diamonds are methanol 
masers (associated wi th star-forming regions ) take n from iGreen et alj d201 lb . Asterisks represent Hn regions observed in radio 
(lAnderson et al Cluster members from lD07l are shown as empty circles. 
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Fig. 9. Run of the ext i nction in several directions covered by our survey. The data have been obtained by applying the technique of 
Cabr era-Lavers et ail J2005) to the red clump giant population within 15' of the field centre. The left panel shows the results for the 
field centred on Stephenson 2. The middle panel represents a field centred on (I = 26°5, b = -(f.5). The right panel corresponds to a 
field centred on (I = 27?5, b = -(f.5). 



Its interpretation is difficult. The dynamical distance for this ve- 
locity range is ^4 kpc, but the extinction does not seem to in- 
crease significantly at this distance. We could think that this 

G26.644+0.020b (v LSR =72.9 kms 4 ), G26.721+0.170b 
(vlsr=70 km s- 1 ) and G26.863-0.275a (v LS r=69.3 km s" 1 ). 



structure is tracing the far side of the Scutum-Crux arm. Other 
possible interpretation would assign (some of) the H n regions 
to the 3 kpc arm, though the RSGs in this velocity range seem 
to have too low reddening to be that far. Finally, we could also 
speculate with the possibility that this feature is due to a strong 
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systematic deviation from the Galactic rotation curve in this di- 
rection. 

Finally, a few H n regions[3 in the neighbourhood of Ste 2 
present velocities around vlsr = +90 kms~'. There are several 
supergiants with vlsr close to these values[3- All these objects 
have E(J - K$) ^ 1.0, suggestive of high distances. 

There is thus some evidence for a small difference between 
the sightline to Ste 2 and neighbouring ones, with a clear lack of 
spiral tracers in the vlsr~ 40 - 60 kms -1 range, in sharp con- 
trast to the I = 24° - 26° range, where H n regions and molec- 
ular clouds are very numerous, and the I- 28° - 30° range, 
where they are found in moderate numbers. The relative distri- 
butions of these tracers may explain the very different redden- 
ings to the RSG clusters in the Scutum Complex. Conversely, 
just in front of Ste 2, there is a fair number of tracers around 
vlsr ~ 70 kms" 1 , not found in the neighbouring regions (see 
Fig.l8land lRathborne et al.ll2009l) . 

5.2. Nature and extent of Stephenson 2 



|D07| claimed 26 likely members of Ste 2 by selecting stars with 
100km s~' < vlsr < 120kms _1 . This choice resulted in an aver- 
age radial velocity for the cluster of vlsr — + 109.3 ±0.7 km s~\ 
The velocity range chosen is somewhat arbitrary. Indeed, D23, 
with vlsr = +119.3 kms -1 is a clear outlier. Star D17, with 
vlsr = +101.4 km s is also far away from the cluster average, 
and has a CO bandhead EW well below that of all other mem- 
bers. The other 24 objects, however, form a very compact group 
in velocity space with 104 kms -1 < vlsr < 113kms _1 . 

Is this ran g e co mpatible with a single cluster? 
iMermilliod et alJ |2008) measured accurate radial veloci- 
ties for several RSGs in Northern clusters, finding aperiodic 
variability in all of them, with variable amplitude in function 
of time. The full amplitudes observed ranged from ^ 5 to 
10 kms 1 . In view of this result, it is fully expectable to find 
a given RSG with v ra d up to 5 kms -1 away from the cluster 
average, even if we do not consider the cluster's dynamical 
velocity dispersion. M oreover, all the RSGs surveyed by 
Mermilli odetail {2008) have moderate luminosities and rather 
early spectral types. It is interesting to note that the two stars 
with late st spectral types, V589 Cas (M3 lab-lb) and RS Per 
(M4 lab: Humphreys!! 19701) . display the largest amplitudes. It is 
possible to envisage later-type M supergiants, many of which 
are dust-enshrouded, as the high-mass analogues of extreme 
AGB stars, subject to heavy pulsation. Objects like VX Sgr 
(M4e -M10ela) o r S Per (M3Iae) display changes in spectral 
type (Wing 2009) and irregular photometric v ariations with ful l 
amplitudes in excess of 5 mag in the V band dKiss et al . 2006). 
Therefore even larger variations in v ra d are n °t discarded. 

Two very bright stars associated with IRAS sources, Dl and 
D4, were not considered members by DO70. a s they both had 
vlsr ~ +95 kms 1 . On the other hand. Deguc hi et alJ d2010al) 
observed masers associated to five sources in the area. For 
D2, they measured velocities vlsr~ +102.7 km s _I and vlsr~ 
+ 106.8 kms 1 in the SiO / = 1-0 v = 1 and H 2 6 16— 523 transi- 



14 We find G026.433+00.614 (vlsr = 89±5km s -1 ), G026.536+00.416 
(v LSR = 87 ± 5 km s-') and G026.261+0.280a (v LSR =88.6 ± 0.4 km s" 1 ). 

15 For instance, J18395861-0549549 (G5Ia), J18402277-0521377 
(M0 la), J18421794-0513098 (M7I-II), J18432089-0431281 (K3Iab) 
or J18433896-0510524 (M7I-II) 

16 Though D07 argue that, given its extreme IR luminosity, Dl could 
be a cluster member experimenting some sort of expansion phase asso- 
ciated to extreme mass loss. 
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Fig. 10. Spatial distribution of candidate red luminous stars, 
passing all the photometric criteria listed in Section [2] in the 
immediate vicinity of Ste 2. Blue circles represent stars with 
i < 15.1, while red squares represent candidates with i > 15.1 
or not detected in i (and so too faint to be observed with AF2). 
Filled circles mark stars actually observed. Members of Ste 2 
from are marked by crosses (note that a few of them have 
not been selected by our criteria). 



tions, respectively. For D49, they measured vlsr~ + 102.7km s _1 
in the SiO transition. Finally, for Dl, they measured vlsr~ 
+92.7 kms -1 and vlsr= + 95.6 kms 1 . The small variations 
with respect to the data in |D07| are likely attributable to the 
very different layers being measured. Degu chi et al.l d2010al) ar- 
gue that the group of sources to the Southwest of the cluster 
(Dl, D5 and D49, plus another source outside the area sur- 
veyed bv lD07l which they call St2-26, and perhaps some others) 
are too far away from the cluster core to be physically related 
and must form a different cluster, which they name Ste 2 SW, 
with some (but non-significant) evidence for lower vlsr- They 
also find vlsr= +86.7 kms 1 in the H2O 616— 523 transition 
(their only detection) of D9. This is very different from the 
vlsr= +111.8 km s _1 measured bv lD07l 

Our data show that the distribution of RSGs in the vicinity of 
Ste 2 is very elongated in the NE-SW direction, forming a nar- 
row structure at least 15' long. Is this very strange shape due to 
an observational bias? This is unlikely, as the DSS2 images sug- 
gest the presence of dark clouds precisely in front of the cluster, 
with significantly lower stellar density than a few arcmin to the 
East. To investigate this issue, in Fig. [10] we plot all the stars 
that pass the colour cuts in 2MASS in the immediate vicinity 
of Ste 2, identifying those that are sufficiently bright in / to be 
observable in our survey. The distribution of sources fully rules 
out the possibility that the known membership of Ste 2 is deter- 
mined by a selection bias. The cluster is surrounded by an area 
of low density of candidates in all directions except the North 
East. In addition, none of the stars observed immediately to the 
West of the cluster has spectral types or v ra d similar to cluster 
members (e.g., D12, D38, J18384958-0557222), while most of 
the candidates observed to the NE are red supergiants with v ra d 
compatible with an association. 

Further to the Northeast, we find the area with a high 
conce ntration of RSGs studied by iLopez-Corred oira et al. 
(11999b . centred close to (/= 27°0, b= +0°1). Though 
these authors do not report accurate coordinates for their ob- 
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jects, we may have some targets in common. For example, 
TMGS J184217.8-051306 (which they classify M8 1) is almost 
certainly 2MASS J 1842 1794-05 13098 (M7I-II here). Pending 
the determination of radial velocities for RSGs in this area, it 
is sensible to assume that at least a fraction of them will have 
similar colours and reddenings to Ste 2 members. 

The Hn region closest in the sky to Ste 2 is 
G026.103-00.069b, just to the South of the cluste r. This 
region has vlsr = 104 + 2 kms 1 d Anderson et al.ll2~01 lb . There 
are several Hn regions in the neighbourhood with similar 
velocities, for instance, G025.766+00.21 1 G025. 800+00.240 or 
G026. 555-00. 305. In principle, these velocities are compatible 
with that of the cluster. Small offsets between tracers probing 
different components are expected, as tracers of H n generally 
see the velocity of the expanding gas. Comparison of velocities 
derived from Ha, a. tracer of the ionised gas, to velocities 
derived from CO, a tracer of the cold gas, for a large sample 
of Hn regions and associated molecular clouds gives a mean 
offset as 4 kms~' (in the sens e of higher CO v elocities), with 
a dispersion of « 12 km s" 1 (jFich et alj |1990(). On the other 
hand, using radio spectra, I Anderson & Bania ([2009) place all 
these H n regions in the vicinity of Ste 2 beyond corrotation, at 
a distance ~ 9 kpc. 

Is there a physical relationship between all these objects 
then? The radial velocity of Ste 2 is compatible with the terminal 
velocity in this direction. Even though some Galactic radial ve- 
locity curves predict higher terminal velocities, tracers that can 
be seen through very heavy extinction, such as H n regions and 
methanol masers, do not reach higher velocities. In view of this, 
we cannot guarantee that every single tracer with velocity around 
vt. sr = +110kms~ 1 is spatia lly close to Ste 2. Indeed, according 
to lAnderson & Banial d2009l) . the Hn regions in its immediate 
neighbourhood with similar vlsr are ~3 kpc away. 

The RSGs observed in this survey, however, are affected by 
extinction. The extinction calculated for several lines of sight 
(Fig. [9j shows a phenomenal rise at d S 6 kpc, reaching values 
E(J -K$) > 2.5. Our sample only contains one star with compa- 
rable reddening, J 18424927-0445095 with E(J - K s ) = 2.8. 
This is a source from the control region with / = 16.1. The 
vast majority of sources with radial velocities similar to Ste 2 
have E(J - Ks) < 2, compared to a cluster median value 
E(J - K$) = 1.76. Therefore none of these RSGs may be be- 
hind the cluster by any substantial amount. On the other hand, 
their LSR velocities force them to be located at distances not 
much shorter than ~ 6 kpc. In view of this, we have to conclude 
that, even though a physical relation (i.e., membership in a single 
association) cannot be strictly claimed with the data available, 
all the RSGs that we detect with vlsr ~ +110 kms -1 must be 
concentrated in a small region, unlikely to be more than 500 pc 
deep. 

Further analysis of this issue is complicated by the strong 
dependence on two basic assumptions: 

i) Kinematic distances assume no strong deviations from the 
Galactic rotation curve. Some of our results suggest that this 
assumption may break down along this line of sight. 

ii) Other distance estimates to Ste 2 assum e an extinction law 
in the infrared close to the standard law dRieke & Lebofskvl 
Il985h . Even though there have been suggestions that this 
law is almost universal (Indeb etouw et al.l 120051) . recent 
work has convincingly shown that the extinction law along 
high-reddening directions towards the inn er Galaxy implies 
much higher transparency in the band |Nishivama et alj 
120081: [Stead & Hoard 120091: iGonzalez et al.ll2012l) . If this 



new law was to be adopted, the distance to Ste 2 would 
certainly increase. 

Testing these two assumptions is clearly outside the scope of 
this paper. Meanwhile, considering the data currently available, 
we can conclude that Ste 2 consists of a compact cluster with ra- 
dius between 2' and 2'.6 - 3.2 to 4.2 pc at 5.5 kpc -, which would 
contain between 16 and 20 RSGs (depending on the radius cho- 
sen), and a much larger halo of stars with similar parameters, a 
large fraction of which are very likely part of a physical associ- 
ation. Even though we cannot be conclusive about the parame- 
ters of the extended association, it likely contains several dozen 
RSGs. Accord ing to current estim ates based on population syn- 
thesis models (Clar k et all200 9b). the central cluster would con- 
tain between ~ 3 and 6 x 10 4 Mq, and so be gravitationally 
bound. The association surrounding it would be much more mas- 
sive, and likely contain many small clusters, which are difficult 
to identify because of the large number of unrelated RSGs pro- 
jected along the line of sight. 

Could the whole overdensity of RSGs be a single associ- 
ation? To answer this question, we need first to determine its 
actual boundaries. With the current data, this possibility cannot 
be ruled out. The most distant RSGs in our sample, lying ~90' 
away from Ste 2, would be located 145 pc away at a distance of 
5.5 kpc. Such size is comparable to that of t he largest star form- 
ing co mplexes in the Galaxy, such as W43 {Nguyen Luong et all 
1201 ll) and W51 dNanda Kumar et alJ 120041) Such complexes, 
however, do not result in a single starburst of star formation, but 
rather proceed to sequential (but not necessarily trig gered) star 
form ation at multiple sites throughout the cloud (e.g., Clark et al. 
2009b). Their end products may be larg e associations, such as 
the complexes of open clusters in M51 (Bastian et al. 2005b or 
the 30 Dor extend ed association, surro unding the central star- 
burst cluster R136 (Wal born et al.ll2002l) . which display a wide 
age range (~ 5 Myr). Conversely, a moderately massive associ- 
ation, such as Per OBI, with a mass « 2 x 10 4 Mq in the core 
region, doe s not seem to show strong evidence for a significant 
age spread dCurrie et al .l l2oToh . 

With an age ^ 20 Myr, the Ste 2 association has had the time 
to evolve dynamically. Therefore the actual distribution of RSGs 
could reflect both a widely multi-seeded star formation process, 
with many progenitors of RSGs forming in small clusters, or the 
consequences of dynamical evolution, with many progenitors of 
RSGs being ejected from Ste 2 and other clusters. If Gaia could 
provide us with accurate proper motions for a sizable sample 
of our RSGs, it could be possible to reconstruct the dynamical 
evolution of the complex and favour one of the two options. 

6. Conclusions 

In this paper, we have presented the results of a spectroscopic 
survey of red luminous stars along the Galactic Plane in the 
/ = 26°-28° region, covering the vicinity of the starburst cluster 
Stephenson 2 and a control region V. 1 away. We have observed 
~250 stars amongst > 600 candidates (> 400 actually observable 
with our instrumentation), finding that all except one are red lu- 
minous stars. We find ~80 red supergiants, most of M spectral 
type. This number represents a significant fraction of the RSGs 
known in the Galaxy, in spite of the fact that we have only sur- 
veyed ~1.7 deg 2 , and checked ~40% of the candidates. The ca- 
pability to select red luminous stars with very simple colour cuts 
has been fully demonstrated. 

The RSGs found are not strongly concentrated towards Ste 2 
(we stress that our observation technique does not allow us to 
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sample but a very small fraction of densely concentrated pop- 
ulations), but seem to be rather uniformly spread over the area 
surveyed. Moreover, about half the RSGs have radial velocities 
and reddenings that make a physical connection with the cluster 
very unlikely. 

On the one hand, these findings demonstrate beyond 
any reasonable doubt that the clusters of RSGs are not 
isolated, but rather part of extended stellar associati ons, a 
result that had already been found f or RSGC3 dPaper It 
iGonzalez-Fern andez & Negueruela! 120121) . This diffuse popu- 
lation likely extends over several squared degrees and com- 
prises several hundred objects. Based on their reddenings and 
radial velocities, members of this population cannot be lo- 
cated at distances from the Earth significantly different from 
those of the RSGCs. Perhaps they represent the result of star- 
burst caused by the tip of th e Galactic bar, as suggested by 
lLopez-Corredoira et al.l (Il999h . Wider-field searches can con- 
firm the existence of such a large population and probe its spatial 
extent. 

On the other hand, the high number of RSGs not associ- 
ated with Ste 2 bears witness to the richness of the young stellar 
population towards the inner Galaxy. Even though the tip of the 
Galactic bar may create special conditions at the distance of the 
Scutum Complex, the foreground population along this sightline 
is unlikely to be affected. If so, such high numbers could be rep- 
resentative of all sightlines towards the inner Milky Way. 

Our results show that the limits of Stephenson 2 cannot be 
clearly defined. A compact core, with radius 5= 2.'5 and contain- 
ing about 20 RSGs, is surrounded by an extended association 
that merges into a general background with an important over- 
density of RSGs. To investigate the spatial extent of this over- 
density and test whether it covers the whole region containing 
the clusters of red supergiants (I ~ 24° -29°), a wider-field spec- 
troscopic survey is needed. Such a survey is currently underway 
and will be presented in a future paper. 

Finally, the technique used here offers a novel way of trac- 
ing Galactic structure. Red supergiants are very bright objects 
in the infrared and can be seen through heavy extinction. Even 
though some of them may have peculiar velocities, their large 
numbers and their membership in clusters make them very good 
tracers. Small clusters containing ~5 RSGs must be relatively 
common in areas of high young population, and they can serve 
as very good tracers of Galactic structure. Even if they are too 
reddened to be accessible around the Ca n triplet, they are very 
good targets for moderately-high-resolution observations in the 
K band. Investigation of further photometric criteria, making use 
of a combination of 2MASS and GLIMPSE colours, seems the 
next natural step along this approach to penetrating into the inner 
Milky Way. 
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